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We show that the combination of a fluid with a bulk dissipative pressure and quintessence matter 
can simultaneously drive an accelerated expansion phase and solve the coincidence problem of our 
current Universe. We then study some scenarios compatible with the observed cosmic acceleration. 



I. INTRODUCTION 



In the wake of the recent measurements of distant exploding stars, supernovae, the existence of negative-pressure 
ON ' dark energy has begun to gain broad consideration. Using Type la supernovae as standard candles to gauge the 
expansion of the Universe, observers have found evidence that the Universe is accelerating JjJ-^]. A new component 
' with significant negative pressure, called quintessence matter (Q-matter for short) will in fact cause the cosmic 
expansion to speed up, so the supernovae observations provide empirical support for a new form of energy with strong 
negative pressure 

, Different forms for the quintessence energy have been proposed. They include a cosmological constant (or more 
■ generaly a variable cosmological term), a scalar field 0, a frustrated network of non-Abelian cosmic strings, and a 
frustrated network of domain walls jl0| , jllj . All these proposals assume the Q-matter behaves as a perfect fluid with 
a linear baryotropic equation of state, and so some effort has been invested in determinating its adiabatic index at the 
present epoch -see e.g. [p|Jl^]. This new energy is to be added to the more familiar components: i.e., normal matter 
(luminous and dark) plus radiation. The contribution of the radiation component is known to be negligible at the 
' current epoch whereas the main contribution to the former comes from cold dark matter (CDM) (b|. 

However, QCDM models (including those in which the quintessence energy is just the energy of the quantum 
vacuum) find dificultics in explaining why the energy densities of the CDM and Q-matter should be comparable 
today. Since both energies redshift at different rates the conditions at the early universe have to be set very carefully 
c/5 for both energy densities to be of the same order today, though one may always invoke some version or other of the 
anthropic principle to soften a bit the problem. This is the coincidence problem Jl4|. Recently, a promising solution, 
for spatially-flat metrics, based in the notion of "tracker fields", fields that roll down their potential according to an 
attractorlike solution to the equations of motion has been proposed. Unfortunately it is not clear to what 

extent these fields have the ability to solve the coincidence problem and, at the same time, drive the Universe to the 
current phase of accelerated expansion. 

In any case, all these models overlook the fact that since the cosmic fluid consists in a mixture of different fluids 
a dissipative pressure may naturally arise which, for expanding universes, is bound to further decrease the total 
pressure Recently it has been proposed that the CDM must self-interact to explain the structure of the halos 
of the galaxies jL7j -see however |l8| . This self-interaction leads naturally to a viscous pressure whose magnitud will 
depend on the mean free path of the CDM particles. On the other hand it has been suggested that dissipative fluids 
(or equivalently particle production processes) can drive a phase of accelerated expansion -see fl9f| and references 
therein. 

This paper investigates how the combined action of dissipative normal matter and a quintessence scalar field may 
lead to the current accelerated expansion stage, and at the same time provide a solution to the coincidence problem 
different from that relying on a tracker field. It is shown in Sec. II that the flat coincident solution is an attractor. 
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Section III deepens in the consequences brought about by a dissipative pressure in the strees-energy tensor, i.e., it 
explores the dynamics of quintessence-dissipative dark matter (QDDM) models, while Sec. IV briefly outlines some 
specific models. Finally, Sec. V summarizes our findings. It should be understood that this work mainly deals with the 
present and late Universe since it is precisely at these stages where the dynamic effects of quintessence and dissipative 
pressure become important. Units in which c = 8ttG = fcs = 1 are used throughout. 



II. COSMOLOGICAL PROBLEMS 

This section shows that the Friedmann-Lamaitre-Robertson- Walker (FLRW) universe filled with perfect normal 
matter plus quintessence fluid, corresponding to some scalar field governed by Klein-Gordon equation, cannot at the 
same time drive an accelerated expansion and solve the coincidence problem. To solve it, without abandoning the 
FLRW geometry, some additional contribution to the stress-energy tensor, such as a bulk dissipative pressure, is 
needed. 

The overall stress-energy tensor of the cosmic fluid without the dissipative pressure reads 

Tab = PU a U b + ph ab , (h a b = 9ab + U a U b V, a U a = -l), (1) 

where p = p m + and p = p m + p^. Here p m and p m are the energy density and pressure of the matter whose 
equation of state is p m = (7™ — l)p m with adiabatic index in the interval 1 < 7™ < 2. Likewise p^ and p^, the energy 
density and pressure of the minimally coupled self-interacting Q-matter field 4>, i.e., 

= \tf + V(<j>) , P4> = \fi - V(<t>) , (2) 

are related by an equation of state similar to that of the matter, viz. p§ — (7^ — l)/9</>, so that its adiabatic index is 
given by 

where for non-negative potentials V(4>) one has < 7^ < 2. The scalar field can be properly interpreted as Q-matter 
provided 7^ < 1 -see e.g. As usual an overdot means derivative with respect to cosmic time. In general 7^ varies 
as the Universe expands, and the same is true for 7 m since the massive and massless components of the matter fluid 
redshift at different rates. 

The Friedmann equation together with the energy conservation of the normal matter fluid and quintessence (Klein- 
Gordon equation) are 

k 1 
3 ( 



H 2 + — = -(p m +Pcb) {k = 1,0,-1), (4) 



pm + 3Hj m p m = 0, (5) 

4> + 3H<j) + V' = 0, (6) 

where H = a/a denotes the Hubble factor and the prime means derivative with respect to <j>- Introducing Q m = p m / p c , 
0^,= p^/pcj with p c = 3-ff 2 the critical density, and flk = —k/(aH) 2 plus the definition fl = Q m + ft^ the set of 
equations can be recast as (cf. |po| ) 

Cl m + fi + fi fe = 1, (7) 

fi = fi(fi-l)(37-2)il, (8) 

fi = [2 + (37 - 2) fi - 3 70 ] Q^H , (9) 

where 7 is the average adiabatic index given by 

7fi = 7 m fi m + 70fi0. (10) 

Next subsections investigate the flatness and coincidence problems. The former will be solved if the solution fi = 1 
to equation (g) becomes an attractor at late time. In its turn, the coincidence problem will be solved if the ratio 
fi^/fim becomes asymptotically a constant. We shall therefore explore the possibility of constant stable solutions in 
the (fi, fi m ,fi^) space. 
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A. The flatness problem and accelerated expansion 



The combined measurements of the cosmic microwave background temperature fluctuations and the distribution 
of galaxies on large scales seem to imply that the Universe may be flat or nearly flat |Q,|l|,^2) . Hence the interesting 
solution at late times of (||) is £1 = 1 (i.e., k = 0), and so we discard the solution £1 = as incompatible with observation. 
The solution £1 = 1 is asymptotically stable for expanding universes (H > 0) provided that the condition dti/dfl < 
holds in a neighborhood of CI = 1 and this implies 7 < 2/3. Hence the matter stress violates the strong energy 
condition (SEC) p + 3p > and as a consequence the Universe accelerates its expansion, i.e., a/a = — (p + 3p)/6 > 0. 

Let us examine more closely the implications of the current accelerated expansion for the QCDM model. Since the 
mixture of Q-matter and perfect dark matter fluid violates the SEC, 7^ must be low enough. Namely, since 7 < 2/3, 
7m > 1, and 7^ < j m , equation (|lo|) implies 7^ < 7. Then, introducing £1 = 1 in equation (^) we obtain 

f} = 3(7 - -y^H, (11) 

and therefore 0^ > 0, i.e., SI4 will grow until the constraint (Q) is saturated, giving = 1 in the asymptotic regime. 
Thus the matter fluid yields a vanishing contribution to the energy density of the Universe at large time. This 
implies that a flat FLRW universe driven by a mixture of normal perfect fluid and quintessence matter cannot both 
drive an accelerated expansion and solve the coincidence problem. Therefore some other contribution must enter the 
stress-energy tensor of the cosmic fluid, i.e., it must be "enlarged". 



B. The coincidence problem 

As shown above we cannot have both the current accelerated expansion and the coincidence problem solved within 
a model that assumes a perfect matter fluid and a quintessential scalar field. Things fare differently when a dissipative 
pressure enters the play. Because of the FLRW metric, velocity gradients causing shear viscosity and temperature 
gradients leading to heat transport must be absent. Therefore the only admissible dissipative term corresponds to a 
bulk dissipative pressure tt. This quantity is always negative for expanding fluids (i.e., n < so long as H > 0) and 
may be understood either as a viscous pressure or as the effect of particle production. On general grounds the former 
possibility is usually thought to give just a small contribution to the overall pressure, however the impact of the latter 
is not so much limited. The expression of the bulk stress when interpreted that way is n = — (p m + p m )T/3H where 
r denotes the particle production rate. This process is dissipative in the sense that the produced particles imply an 
augment of the phase space volume. A recent discussion about the interplay between dissipative bulk pressure and 
cosmological particle production can be found in l|J . 

So the total stress-energy tensor of the cosmic medium, made up of a dissipative but otherwise normal fluid plus 
the Q-matter fluid, reads 

+ p<p + ir)u a u b + (p m +Ptf, + ir)g a b- (12) 
A parallel calculation to that of above leads to the corresponding Einstein-Klein-Gordon field equations 



n = fi(o- 1) 



3( 7 +-)-2 



(13) 



and 



3(7+- I -2 



£1 - 370 > HCl^, 



(14) 



instead of equations fla) and (Bl) . The energy conservation of the normal matter is 



Pm + 3 ( 7m H ) PrnH = 0. 

Pn 



(15) 



Owing to the presence of the dissipative bulk stress the constraint 7 < 2/3 does not longer have to be fulfilled for 
the solution £7 = 1 of equation ( |l3| ) to be stable. Likewise, inspection of (|lj) shows that when £1=1 one can have 
tl^ < just by choosing the ratio tt/ p sufficiently negative . Thereby the constraint (0) allows a nonvanishing £l m at 
large times. By contrast tracker fields based models (valid only when £lfe = 0) predict that fi m — > asymptotically 
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P- 5] . A fixed point solution of equation (y_3J) is f2 = 1. Note that equations (m) and (|14|) have fixed point solutions 
Q m = £l m o and = f^Oj respectively, when the partial adiabatic indices and the dissipative pressure are related by 



TT 

7m H 

Pm 



14> 



2H 
3IP 



(16) 



Then smaller 70, the larger the dissipative effects. Let us investigate the requirements imposed by the stability of 
these solutions. From ( |l3| ) we see that 7 + tt / p < 2/3 must be fulfilled if the solution £1 = 1 is to be asymptotically 
stable. This condition, together with ([fq), leads to the additional constraint on the viscosity pressure 



TT < 



m I ym , 



(17) 



which is negative for ordinary matter fluids. Also by virtue of ([Tl]) and the first equality in we obtain from ( |17j ) 
that 70 < 2/3. 

In the special case of a spatially flat universe (f2 = 1), the stability of the solutions fimo and fi^o may be studied 
directly from (Q). Namely, setting = f^o + oj and using ([[o]) it follows that 

(18) 



Co = 3f2 m I 7™ - 70 H I H (Q$o + w) . 

V Pm J 

Accordingly the solution fi = 1, fi^ = O^q is stable for the class of models that satisfies ip = j m — 70 + (n/ p m ) < 
and ip — > for t — > 00. Note that this coincides with the attractor condition (|l6|). 

In order to study the stability of the solutions Q m o and $7,00 when k ^ it is advisable to derive a dynamic equation 
for the density ratio parameter 

^ 

^0 

For this purpose we combine the logarithmic derivative of e with the definitions of f£ m and Q$ and the energy 
conservation equations ( |l5| ) and (|^) -the latter written in terms of p^. It yields 



3 70 - 7r> 



To calculate 70 we use (|15[) together with Q3J), (Ml) and (H), obtaining 



70 = 7m + 



2# 7T 

3fP 



Introducing (EOl) in (19) we get 



3iJe 



2i? TT 
—72 + 7 m + — 



3ff 2 



o - 7m S2fe 

3 p, 



(19) 



(20) 



(21) 



and perturbating this expression about the solution eo 
obtain with the help of ([[ 
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O(l), (i.e., using the ansatz e = eq + 5 with \S\ <C 1) we 



(22) 



70 Q k H (e + 5) 



near the attractor. For flf. > (negatively spatially curved universes) it follows that 6 decreases, i.e., the ratio 
(f2m/^0)o is a stable solution. For f2^ < one has to go beyond the linear perturbative regime and/or restrict the 
class of models as in the spatially flat case to determine the stability of the solution. We defer this to a future research. 

As we mentioned above, recently there have been some claims that CDM must not be a perfect fluid because it 
ought to self-interact (with a mean free path in the range 1 kpc < I < 1 Mpc) if one wish to explain the structure 
of the halos of galaxies |17| . In this light it is not unreasonable to think that this same interaction is the origin of 
the dissipative pressure tt at cosmological scales. Bearing in mind that I — 1/na, with n the number density of CDM 
particles and a the interaction cross section, a simple estimation reveals that at such scales I is lower than the Hubble 
distance H^ 1 and accordingly the fluid approximation we are using is valid. 
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III. QDDM ASYMPTOTIC ERA 



Bulk viscosity arises typically in mixtures-either of different particles species, as in a radiative fluid, or of the same 
species but with different energies, as in a Maxwell-Boltzmann gas. Physically, we can think of ir as the internal 
"friction" that sets in due to the different cooling rates in the expanding mixture. 

Any dissipation in exact FLRW universes have to be scalar in nature, and in principle it may be modelled as a 
bulk viscosity effect within a nonequilibrium thermodynamic theory such as the Israel-Stewart's |24| , |2j|. In that 
formulation, the transport equation for the bulk viscous pressure takes the form 



TTT 



1 



3H 



T 
T 



(23) 



where the positive-definite quantity £ stands for the phenomenological coefficient of bulk viscosity, T the temperature 
of the fluid, and r the relaxation time associated to the dissipative pressure -i.e., the time the system would take to 
reach the thermodynamic equilibrium state if the velocity divergence were suddenly turned off |26| . Usually C is given 
by the kinetic theory of gases or a fluctuation-dissipation theorem or both [^7| . 

Provided the factor within the square bracket in ( ^3|) is small it can be approximated by the more manageable 
truncated transport equation 



-3CH, 



(24) 



widely used in the literature. This as well as (|23|) meets the requirements of causality and stability to be fulfilled by 
any physically acceptable transport equation |2gjj. 



A. The quasiperfect regime 



Here we obtain an explicit expression for the leading behavior of the attractor solution at late time. We begin by 
writing the equation of motion for the Hubble factor that follows from combining ([l^) and ( |l5| ) with ( pi| ) 



H + 3~fHH + t 
k 



[(l-§ 7 ) (2i/-r- 1 ) + | 7 ]=0. 



(25) 



We next evaluate (10) on the attractor and insert it together with (|16[) in (J23) to obtain 




37m H 



H 



3 7 „ 



■H 2 



3C 



20 



H = 0. 



(26) 



Observation seems to rule out huge entropy production processes on large scales, otherwise the flux of gamma-rays 
we witness should be much higher j29|. Hence we shall assume that the viscous effects are not as large as that, but 
however not altogether negligible. If r is the relaxation time, then v = (tH) -1 is the number of relaxation times in 
a Hubble time -for quasistatic expansions v is proportional to the number of particle interactions in a Hubble time. 
Perfect fluid behavior occurs in the limit v — > oo, and a consistent hydrodynamical description of the fluids requires 
v > 1. Thus we are lead to assume that tH is small and we propose a "quasiperfect" expansion in powers of v~ l . 

Let us show that the attractor solution of leading behavior a ~ t a when t — > oo, with a a positive-definite constant, 
is consistent in the quasiperfect regime. Indeed, by virtue of (|l6|) it implies 7^ ~ 2/3ct and 



For approximately constant 7m ( 7r , 



7T / 2 

— ~ - 7 m - — 

Pra \ 3<7 

= 1 for CDM), we get from (| 



(27) 



Hence (24) becomes 



7T 1 



2 

va 



and we get to leading order in v 1 



C - ^mO ( 7m - — ) H . 



(29) 



(30) 



Then integration of ( |l5| ) yields p m ~ a -2 /' 7 , the same scaling law as Finally its insertion in (||) leads back to 
a ~ i" 7 , showing the consistency of our assumptions. These results correspond to the lowest order in the quasiperfect 
expansion. 

To go a step further we introduce the expansion of H in powers of v^ 1 



H = H (l + h x v- x 



in (26), and assuming that \f\ -C v , it follows the approximated solution 



t 



1 



37 m CT - 2 



1 

i / v 



(31) 



(32) 



where is an arbitrary integration constant. This expression reveals that the power law is an attractor solution 
and that for a > 1 (deceleration parameter q = — (a — I) /a < 0), CDM viscosity provides an accelerated expansion 
scenario that also solves the coincidence problem. It can be shown that y$ does not pick any correction of order t~ 1 



from the subdominant term in (|32|). Instead the first correction appears to the order v~ t~~ , and this fact shows the 
high degree of correction of the approximation that 7^ takes a constant value in the late time regime. 

We note that this attractor solution works for any viscosity coefficient with leading behavior ([BO]) . In particular the 
case C oc -Jp m , investigated in |^,^l[, satisfies this requirement. 



B. Full causal corrections 



Here we gauge the changes bring about by the the full transport equation ( p3| ) on the expansion exponent obtained 
in the previous section. Using that equation and the viscosity coefficient found in (p0|) we get 



l + 2rH 2 
2 IP 



7m 1 - - r 



H + - lm H J 



H + -H 2 = 0, 



(33) 



where, to estimate the corrections, we have assumed that in the asymptotic regime T tx p r , with r a positive-definite 
constant, and we have used that p m ~ eop/ (1 + eo) in this regime. This power law relationship is the simplest way 
to guarantee a positive heat capacity. Usually p, p, T and the particle number density n are equilibrium magnitudes 
related by equations of state of the form p = p(T, n) and p = p(T, n). Further the thermodynamic relation 



dp 

On 



p + p T I dp 
n n \ dT 



(34) 



holds. This directly follows from the requirement that the entropy is a state function [p2| . In the particular case of a 
material fluid with p = p(T) and constant adiabatic index, this relation imposes the constraint r = (7 — 1) /7, so that 
< r < 1/2 for 1 < 7 < 2. Inserting ( |3~l| ) in ([53]), and assuming that |f| <C f -1 , we obtain the approximate solution 



2r + 1 
3cr 



C 

T 



(35) 



Comparison with ( p2[ ) shows that except when a ~ 1, the use of the complete transport equation leads a to slighty 
slower rate of expansion at late time. Also, in this regime, the equilibrium temperature decreases as T ~ t~ 2r ~ a~ 2r l a . 
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C. Late Q-matter dynamics 



In virtue of ( |l0| ) the density parameter ratio can be written in terms of the adiabatic indices 

7m - 7 ' 



(36) 



Since 7 m is approximately constant, if 7 — ► 70 in the asymptotic regime when e — ► €q, then 7^ must also approach a 
constant value. Hence from (||) get the constraint 

VW/j? ~ C (37) 

with C > 1 as 70 < 2/3. Potentials that satisfy this constraint have been investigated in |Q and |54j. Then (^) 
becomes 

* 35 ,0. (38) 



j> 1 + 2C 

An interesting potential that meets this constraint is the exponential potential 

V( ( f>) = V Q eM-M), (39) 
where A and Vq are constants. Now, integrating ( |38] ) and using ( |37| ) we get 



In ^ 1 \+2ht 
2 (2 - 70) 



(40) 



Hence slowly rolls down the exponential potential as oc 1/t when t — > 00. Also we find that C ~ (3er — l)/2 with 
7^ ~ 2/3ct < 2/3 on the attractor, irrespective of Vo and A. 

Perfect fluid QCDM models based on the exponential potential are ruled out by observations 0. However we shall 
demonstrate in the next secion that in the realm of QDDM models the exponential potential yields satisfactory results 
without any fine tuning of the parameters. 

IV. QDDM MODELS 

This section explores the dynamical evolution of a universe filled with a viscous material fluid and a quintessence 
scalar field by resorting to models based on simple relationships for the nonequilibrium quantities. This allows us to 
explore the large dissipative regime where the nonequilibrium pressure has a magnitude comparable with the energy 
density. Recently tracker-field models with inverse power potentials have attracted much interest [l5|j35[| . Here we 
will investigate some QDDM models with exponential potentials (|3^) such that for a wide range of initial conditions 
the scalar field settles into an attractor solution that depends only upon a few nonequilibrium thcrmodynamical 
parameters, addressing the coincidence problem. 

A. Linear dissipative regime 

The linear regime C = otH, with a a constant in the interval < a < 1, arises for instance when the coefficient 
of bulk viscosity takes the form of a radiating fluid. We further assume that the number of interactions of a generic 
CDM particle in a Hubble time is larger than unity so that the hydrodynamic regime is respected. Here we investigate 
models with the limiting behavior 7 — > 2/3 in the asymptotic regime. We begin by inserting the ansatz 

7=l(l+X) (41) 
in equation (|25|). It is immediately seen that the latter splits in two equations, namely 

H + {2 + v)HH + v (l - ^aj H 3 = 0, (42) 



7 



and 



H 2 



2H H 



H 2 



x = o. 



Replacing the solution of (^) in ( [4l| ) it follows, 



2 / , a 2 -" 

3{ 1 + b kT^ 



(43) 



(44) 



where b is an arbitrary integration constant. This expression for 7 will be sensible only if it meets the restriction 
< 7 < 2. Then (f42| ) can be transformed into a linear differential equation of second order whose general solution in 
parametrized form is already known ]3(|, In particular a oc t a , where a is the larg est root of v(\ — 3a/2)er 2 — 
(2 + v)a + 2 = 0, is an asymptotic stable solution in the limit t — > 00. Then, using ( Utj ) together with the attractor 
conditions = 1,7 = 2/3 and (|l6|), we find 



2(1 — fi TO o) 
2 3ri m o7m 



(45) 



Hence the quintessence adiabatic index 7^ = 2/3cr depends solely on the dark matter parameters fl m and 7 TO in the 
asymptotic regime. Moreover, a relationship between the dissipative parameters a and v follows from ( p9| ) and (|lq), 
namely 



a = fi m o 1 



7 '"-37 



(46) 



and the requirement a > implies > 2 or ^ > i/ m »n — ^7^. The same condition arises from the requirement that 
7 — > 2/3 when t — > 00. Equations ( |45| ) and ( p6| ) show that a grows with ;y attaining a max — (37 m — 2)f2 m /3(l — O m ) 
in the limit v — > 00. 

As it follows from (45) there will be accelerated expansion (i.e., a > 1) if f2 m o < 2/37 m < 2/3 and accordingly 
we obtain a family of exact solutions describing a QDDM scenario that solves the coincidence problem regardless of 
the value of the spatial curvature. Figure 1 depicts the dependence of 7^ on fi m when 7 m = 1. To make a rough 
estimate of the cosmological parameters in the late time era we assume that our Universe is currently close to the 
asymptotic attractor regime and use the current observational bounds. After Q the combination of low redshift, 
type la supernovae and COBE measurements determines (for a spatially flat universe) the range Q m ~ 0.3 — 0.4 
and 70 < 0.6. From figure 1 it is seen that our linear dissipative model satisfies comfortably these constraints. For 
Q m — 0.3 we get from ( f45| ) a ~ 1.27. This is fully consistent with current estimations of Ht today ]37]], as they 
provide a lower bound for it in a universe that started only recently a phase of accelerated expansion and approaches 
asymptotically the attractor regime. 



B. Viscous speed regime 

This scenario is somewhat more general than the previous one. It arises when the bulk viscosity coefficient is given 
in terms of the speed of the bulk viscous signal v by p5[ 



C 2 

- = V JmPr, 

T 



(47) 



We further assume r related to H by the same expression as before, only that to simplify the calculations we now 
take v constant. Hence (E3) becomes 



h" + 3 7 // + 3Vh - 9^^h 
1 + e 



+ h' + 3jh = 0, 



(48) 



where 



h = H 2 



(49) 
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and the prime indicates derivative with respect to i] = In a. Here we have used the scale factor a(t) as a coordinate 
instead of the cosmological time t, i.e., a is assumed to be a monotonic function of t. 

Equation is useful to study the asymptotic stabilit y o f FLRW expansions at late time because it can be 
rewritten in terms of the derivative of a Lyapunov function |3q| 



dr)\2 



7' + 71/ — 3 



1 + e 



h 2 \ = - (3 7 + v) h' 2 



3 

+ 2 



7' + — 3 



u 2 7n 



1 



(50) 



If the adiabatic index does not decreases too fast in the attractor era, a sufficient condition for the Lyapunov function 
to have a minimum at the phase space point (h, h') = (0,0) is that v > 3v 2 . Within this scenario the parameters of 
the matter fluid can be taken as quasistatic, yielding an asymptotically stable minimum. The leading behavior of the 
solutions in this quasistatic regime is given by 



cia Al + c 2 a A 



where 



Ai, 2 = -i-(3j+v)± 



(37 - v) + 367 m v 2 fl r 



1/2 



(51) 



(52) 



and the parameters are evaluated at the asymptotic attractor era. For large scale factor, equation (|5^) reduces to 
h 2 ~ c\a Xl when —2 < Ai < 0, and we have once again a power-law accelerated cosmic expansion with a = — 2/Ai. 



Likewise equation (^Tj) reduces to h 2 
to linear evolutions at late time. 
Combining Q and @ we get 



CiO 



for Ai = —2, and h 2 ~ for A x < —2. These two latter cases correspond 



A 2 + (3 7m + v) A + 3 7m (i/ - 3v 2 ) = 



(53) 



and using ( [Tof ) together with the attractor constraints f2 
Hence we obtain 



1 and (O) we find that Ai = —37^ also satisfies (|52 



!<$> 



G 



37n 



(37. 



1/2 



(54) 



In this model the quintessence adiabatic index does depend on the parameters v and v while it is independent of the 
density parameter J7 m . We plot 70 in figure 2 for 7m = 1. As it can be seen a wide range of the parameter space 
(2/, v) is consistent with a spatially-flat accelerated universe and such that 70 < 0.6 for any value of fi m . This shows 
another solution to the coincidence problem for any value of the spatial curvature and compatible with accelerated 
expansion. It is also seen that dissipative effects enlarge the parameter space where observational data has to be 
fitted, but global dynamic information alone cannot determine the specific values of v and v. Note that the smaller 
70, the larger the dissipative contribution to the sound speed, and the smaller the interaction rate. 



V. DISCUSSION 



We have proved that the coincidence problem and an accelerated expansion phase of FLRW cosmologies cannot be 
simultaneously addressed by the combined effect of a perfect fluid and Q-matter. Nonetheless, if nonbaryonic dark 
matter behaves as a dissipative fluid rather than a perfect one, both problems may find a simultaneous solution. This 
is so because an imperfect (i.e., dissipative fluid) expanding in a FLRW background possess a negative pressure 7r 
that enters the conservation equations of general relativity. The models presented here are compatible with a negative 
deceleration parameter at present time. In consequence, the quintessence scenario becomes more robust when the 
dissipative effect of the nonequilibrium pressure arising in the CDM gas is allowed into the picture. 

Recently attempts have been made to constraint the state equation of the cosmic fluids (Q-matter included) by 
considering gravitational lcnsing effects, the mass power spectrum and the anisotropics of the cosmic backgroud radia- 
tion H , f39|| . We have shown specific models with an ample region in the space of out-of-equilibrium thermodynamic 







parameters satisfying this constraint in the asymptotic attractor regime which our Universe may well be approach- 
ing. We would like to point out that the parameter space should be enlarged by adding these out-of-equilibrium 
parameters when fitting the observational data. Unfortunately there is some degenerancy in the determination of 
these parameters from constraints arising from the cosmological dynamics alone. We hope, however, that simulations 
of structure formation that include dissipative effects will ultimately prove instrumental in discriminating between 
different models. 

On the attractor asymptotic regime the dissipative matter fluid and the scalar field contribute in a fixed ratio to 
the pressure and energy density along the QDDM era. This scenario ameliorates the self-adjusting model [[l0| as it 
allows for < 7 < 2 for a wide range of initial conditions. It also improves on the tracking models as it solves the 
coincidence problem in the late accelerated expansion phase. While keeping a finite difference between the adiabatic 
indices of quintessence and matter fluid, this difference arises in the viscous pressure. 

As it has been noted the Q-matter proposal may entail some undesirable effects such as the variation of the 
constant of nature and the presence of unobserved long range forces. Efforts to solve this difficulties by coupling the 
quintessence field to the electromagnetic field pH , and to the curvature of the metric 0] have been made. Further, 
a time dependent but otherwise smooth scalar field such as those studied so far, are somehow unphysical as they 
violate the principle of equivalence -the Q-matter must experience clustering in some degree and so it cannot be 
entirely smooth. Therefore to respect the equivalence principle one should assume that <fi varies with position as well. 
Accordingly one should be led to forsake the FLRW metric an take up some inhomogeneous one instead, only that in 
such a case the computational effort is bound to be enormous and most likely no exact solution will emerge. 

Despite that the Q-matter proposal cannot be regarded at this moment with unreserved confidence, we feel this 
idea is still worth exploring in the hope that the aforesaid difficulties may soon find a satisfactory answer. 
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Figure Captions 

Figure 1. The adiabatic index 7^ of the quintessence scalar field versus the matter density parameter fi m , for CDM 
(7m = 1), in the asymptotic attractor regime for the model presented in Sec. IV A. 

Figure 2. The adiabatic index 7^ of the quintessence scalar field versus the interaction rate parameter v~ x and the dissipative 
contribution to the sound speed v, for CDM (-y m — 1), in the asymptotic attractor regime, for the model presented in Sec. 
IV B. 
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